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Abstract 

We present a new strategy for multipulse control over decoherence. When a two-level system 
interacts with a reservoir characterized by a specific frequency we find that the decoherence is 
effectively suppressed by synchronizing the pulse-train application with the dynamical motion 
of the reservoir. We discuss the applicability of this strategy by studying the dependence of the 
decoherence suppression on the shape of the coupling spectral density. We find that the effectiveness 
of this strategy arises from the non-Mar kovian nature of dynamical motion of the reservoir. 



PACS numbers: 03.65.Yz,03.67.Hk,05.30-d 



I. INTRODUCTION 



Degradation of quantum superposed state by decoherence is an obstacle to quantum 
information processing. In order to proactively prevent errors, a multipulse control method 
has been proposed . It is essential that the application of tt pulses causes time reversal 

in order to suppress the decoherence. The method has attracted considerable attention; 
it has been applied to suppress unwanted spontaneous emission , the magnetic-state 



decoherence by collisions in a vapor 



and the damping of vibrational mode of a chain 



of trapped ions[s[ l|. While the multipulse control method requires no ancillary bits and 
no accurate detection, its effectiveness has been shown when sufficiently short and strong 
pulses are periodical ly appli ed in a shorter interval than the characteristic time of the system- 



reservoir interaction! 10 



ll| . The degree of suppression becomes larger as the pulse interval 



becomes shorter. Since these conditions are not easy to execute, a new approach to use a 
continuous control field instead of pulses has been proposed [12]. It is also shown that the 
control pulses do not always have to be ultra short for systems coupled to the reservoir with 
1/f spectral density jl 31] . In order to overcome the strict condition on the pulse application, 
it is desirable to seek a possibility to use the pulse trains with relatively long pulse interval. 



In previous paper 



14j . we have suggested a possibility to relieve the condition on pulse 



interval by formulating a theory of pulse control on the pure dephasing phenomena that 
is caused by the interaction with a boson reservoir. Since the ordinary spin-boson model 
where a spin linearly interacts with the boson reservoir cannot describe the irreversibility 
in the long time region except for the ohmic dissipation case, we have extended the model 
to include a nonlinear interaction. We have found that the multipulse control is effective 
for this model when the pulse interval is shorter than the reservoir correlation time. We 
also found that the effective pure dephasing time shows a non-monotonic dependence on the 
pulse interval, that is, it has a peak when an application of it pulse-train is synchronized 
with the oscillation of the reservoir. This means that the pure dephasing phenomenon is 
also effectively suppressed by paying attention to the dynamics of the reservoir. However, in 
many cases, the system-reservoir interaction is described with a strong linear interaction and 
a weak nonlinear interaction. In the time region where we want to discuss the effectiveness 
of suppression by synchronizing the pulse application and the dynamics of the reservoir, the 
linear interaction plays an essential role in the decay. 
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In this paper, assuming that a two-level system linearly interacts with a boson reservoir 
that has a characteristic frequency, we discuss the effectiveness of the synchronization of a 
7r pulse train with the oscillation of the reservoir. For convenience, we name this strategy 
as synchronized pulse control(SPC) in the following discussion. As recognized in ^J, SPC 
also depends on the type of coupling (bath) spectral density. In order to make clear the 
applicability of SPC, we study the effectiveness of SPC on non-Lorentzian and Lorentzian 
coupling spectral density. 

The outline of this paper is as follows: In Sec|Tll we introduce the model of the decoherence 
and derive the basic formula for multipulse control on the linear spin-boson model. Next, 
we discuss the synchronized pulse control in Sec lllll The application of the basic formula to 
non-Lorentzian (Lorentzian) coupling spectral density is written in Sec Jllll A.fPA respectively. 
After discussing the effectiveness of the SPC in Sec flVl we give concluding remarks in SecfVl 



II. FORMULATION 

We consider a two-level system composed of an excited state \e) and a ground state \g) 
with energy E e . Let us consider the decoherence of this two-level system, which is caused by 
a linear interaction between the excited state and a boson reservoir. The system Hamiltonian 
reads, 

'Hr = Ho + H-sb = (7~ts + 'Hb) + T~Lsb, (1) 
Us = E e \e)(e\, (2) 
'Hb = h^e k blb k , (3) 

k 

H SB = h\e){e\J2h k e k (b k + bl). (4) 

k 

In order to suppress the decoherence, we apply pulses that are sufficiently short and 
strong. This indicates that the interaction with the reservoir is neglected during pulse 
application: 

N 

H SP (t) = Hs + ^2n P j(t), (5) 

3=0 

H PJ {t) = ■ Jl (\e)(g\e** + \g)(e\e^), (6) 
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where Ej(t) is the j-th applied pulse of external field. We assume the pulse to be on 
resonance with the two-level system, which means E e = Tiuo. 

When we apply N pulses with a pulse interval r s and pulse duration At, the time evolution 
of the density operator pit) of the total system, is given by 

p{t) = e -^fl(*-(^+ A *))T + [ e - i ^ s+Atd *' Lp ^ ( '' ) ] 



N-l 

x ■ - - - - ■■ ■ • r '"" : Ai 



(7) 

where T + is the time ordering symbol from right to left and Lpj (Lr ) indicates the Liouville 
operator during the j-th pulse (the interaction with the reservoir) which is defined as 

iL v ■■■ =\[U V ,•••], (v={P,j}or{R}). (8) 
We rewrite Eq.l jTj l by using the following relation for an arbitrary operator X, 

which consists of the operators as 

e -iH R (t-to) = e -jrH (t-t ) T+ ^ e -y^- to dt'H SB (t')^ 

= (Utit - t )|e><ej + U 2 (t - t )\g)(g\) (10) 

T + [ e -i^o dt ' Hp ^} = e -i H o(t-to) e -iHrAt-to) 

= U 2 (t- t Q ) cos(|) (e— (*"*>) |e> (e| + \g) (g\) 

-i S m( -l)(\e)(g\ + \g)(e\), (11) 
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where 



H SB (t) = e-t Hot H SB e* Hot = ^ h k e k (b k e~^ + b\e^\ (12) 



k 



H P ,it) = e-i a * t H PJ ei B * t = -\e 3 ■ p {\e)(g\ + \g)(e\) = H PJ . (13) 



Here we have assumed the each applied pulse to be square whose strength is Ej , which gives 

M 
h 



pulse area 9j = -j^-(t — t ) for the j-th pulse. In Eq.((TTJ), we used the following definitions 



4 



as, 

l/i(t) = exp[-i(w + 5]e46 fc )*]r + [expH f dt'^hk^he-^' + b\e^ t ')}} 

k Jo k 

= r]{t)e^[-i(u + Y,^blb k )t}eM^ e MUt)-hCk(tm (14) 

C/ 2 (t) = expH^e^M], (15) 

fc 

where 

77(0 = exp[i V>£(e fc t - sine k t)}, &(t) = ^(1 - e ie **). (16) 

Now we suppose the pulse area 9j to be ir except for the first pulse whose pulse area is | 
to generate a superposed two-level state at an initial time (t — 0). Defining the intensity of 
off diagonal element of the density operator p(t) as 

I(t) = \Tr R (e\p{t)\g)\\ (17) 

where Tr R denotes the operation to trace over the reservoir variable, we obtain for even N, 

I(t) = \Tm\U.it - iVr s )(f/ 2 (r s )f/ 1 (r s )) JV / 2 p(0)(f/ 2 t (T- s )f/ 1 t (r s )) Ar / 2 f/ 2 t ( t _ iVr s )]| 2 , (18) 

and for odd N, 

I(t) = \Tr R [U 1 {t-NT s ){U 2 {T s )U 1 {T s ))^ 

(19) 

Here we have assumed the pulse duration At to be infinitely small. In the case of the pulsed 
magnetic resonance or the transient nonlinear optics, I(t) indicates the signal intensity. 

Next, we focus on the time evolution of the boson reservoir, eliminating the two-level 
system that periodically changes its state between |e) and \g) by the 7r pulse train. Denoting 
the displacement operator as 

D({a k })=eMj2(^ b l- a *M (20) 

k 

where {• • • } means a set of bosons in the reservoir, we obtain the off diagonal element of 
the density operator p(t) in the form for even N, 

I(t) = |Tr i? ( e |[(|A JV (t))|e))((5 J v(t)|(^|)]|^)| 2 = \(B N (t)\A N (t))\ 2 , (21) 



5 



with 



\A N (t)) = D({a Ntk (t)})\0) = \{a N , k (t)}) , (B N {t)\ = (0\D({f3 N , k (t)})\ = 

(22) 

Here we defined a^^it) and 0N,k{t) as 

TV N 

a N ^ k (t) = -h k + J2(-^{hke- lek{t - jTs) } , fe(t)E^(-iy- 1 {/i fce ^( t ^) }i (23) 

3=0 3=1 

For odd N, we obtain, 

1(f) = |Tr i? (e|[(|5 J v(t))|e))((A iV (t)|^|)]^)| 2 = \({a N , k (t)}\{p N , k (t)})\ 2 , (24) 

with 

N N 

a ^ k (t) = J2(-^ j {hke- iek{t - jTs) } , P Ntk {t) = ^(-l^-H/ifee-^*-^) - h k }. (25) 

3=0 3=1 

In obtaining Eqs. (j21|) ~ (J2*5j) . we assume that the boson reservoir is in the vacuum state 
and the two-level system is in the ground state at the initial time: 

p(0) = b>fo|®|0><0|. (26) 

Eqs. (121 J) and (}2"l*|) imply that the intensity I(t) is described with overlap between coherent 
states \A N (t)) and \B N (t)) of the reservoir. These coherent states are biuniquely associated 
to \g) and |e): For even N, the state |Ajv(f)) is associated to |e), whereas |-Bjv(t)) is with 
The 7r-pulse application alternately exchange the association between the reservoir states 
and the two-level system. 

An actual evaluation of Eqs.(jU) ~ (|25|) requires us to rewrite the summation over k into 
the energy integral, 

rOO /"OO 

Y, \h k \ 2 fM = V \h k \ 2 f(u k ) / de5(e -u k )= / deh{e)f{e) 1 (27) 
k k Jo Jo 

where we have defined coupling spectral density h(e) as, 

h(e) = J2\ h k\ 2 S(e-^ k ). (28) 
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III. NUMERICAL EVALUATION 



Now we evaluate the time evolution of the intensity I(t) 

I{t) = \(A N (t)\B N (t))\ 2 = exp[-J2 \<*N,k(t) -fe(t)| 2 ]. (29) 

k 

In the following, we discuss the time dependence of I(t) for non-Lorentzian and Lorentzian 
coupling spectral density. While the SPC can be effective for the former case, it is ineffective 
for the latter case. 



A. Non-Lorentzian coupling spectral density 

As the first example of the non-Lorentzian coupling spectral density, we consider a Gaus- 
sian distribution with the mean frequency uo p and the variance j p , 

,iG(e, ^ exp( - l£ ir !) ' (30) 

Setting N = in Eq.(|2I|L we evaluate time evolution of I(t) after a single | pulse at 
t = 0, which is shown in Fig^ Here and henceforth, we have used a scaled time variable 
as t = uj p t and set the parameters as 7 P = 7 p /cu p = 0.15, s = 3, which mean that the decay 
time of the interaction mode is relatively long, the average number of boson which interact 
with the spin is 3. We see a damped oscillation whose period is 27r. 




2n 4n 



FIG. 1: Time evolution of I(t) without pulse control for 7 p = J p /u) p = 0.15, s = 3. 



The dynamical decoupling method [ 



tells us that an application of ir pulse-train 



is sufficiently effective when the pulses are applied in a sufficiently "small" interval. When 
the spectral density has a cutoff frequency at u c as assumed in 3, [ll|, the pulse interval 



t s is required to be much smaller than r c = oo^ 1 in order to control the decoherence. In the 
case of the Gaussian distribution as Eq. fJHUj) with relatively small variance, we suppose that 
the pulse interval is required to be much smaller than io~ x = ^ for decoherence control. In 
FigEJ we can see that the decay of I(t) is well suppressed for r s = 
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1 

FIG. 2: Time evolution of I(t) for the pulse interval r s = Other parameters are same as in 
FiglU 

However, when the pulse interval becomes longer, we find that the pulse application makes 
things even worse than the damped oscillation in FigO This is shown in FigEl where the 
pulse interval is r s = y. Now we plot the case where 7r pulses are applied with the interval 



"1 
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FIG. 3: Time evolution of I(t) for the pulse interval r s = Other parameters are same as in 
Figd 

t s = t p in FigEJ where we find that the phase coherence recovers at the pulse application 
time. The peak value asymptotically goes to be constant, which reflects that the dephasing 
in long time region cannot be described by the linear interaction. When we obtain a recovery 
of the intensity by synchronizing the pulse application with the characteristic period r p , we 
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call this strategy for suppression of decoherence as synchronized pulse control. In this paper, 
we consider only the linear interaction between the spin and the original boson reservoir. 
While it is necessary to take into account the nonlinear interaction, which causes the pur e 
dephasing phenomena (irreversible processes in the long time region) in many systems 1J] , 
the effect of the pure dephasing is not significant in the time region shown in Fig0J since 
the nonlinear interaction is often much weaker than the linear one. 

'1 
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FIG. 4: Time evolution of I(t) for the pulse interval r s = t p . Other parameters are same as in 
Figd 

The synchronized pulse application has been discussed in the context of transient opti- 
cal nonlinear spectroscopy, called synchronized quantum-beat echoes(SQBE) by Tanigawa, 
et.al.jlf]]. They used two light pulse trains in order to cause an optical transition between 
a ground-state sublevel pair. The repetition frequency of the pulse trains is equal to the 
separation of the sublevel pair, which is superficially similar to our approach. However, 
the role of the pulse train in this paper is essentially different from SQBE, because each 
pulse area in SQBE is much smaller than tt. In the pulse train control of decoherence in 
the present work, it is essential that the pulse area of each pulse is tt, because the physical 
origin of the coherence recovery is the time reversal operation caused by each tt pulse. In 
Tanigawa' s work, the maximum signal is generated when the total pulse area of the second 
pulse train is tt, and the pulse train is used to achieve the sublevel resonance. 

Next, we assume the coupling spectral density to have semi-elliptic distribution, 

h s (e) = s^yj-(e - u p ) 2 +p, (31) 

defining p as 

H 



to have half width 7 P . The coupling function has been used to describe the coupling strength 
between phonons and a localized electron in a so lid [l^. 

We show the time evolution of I(t) for the same parameters in Fig0] Fig|5] shows similar 
behavior as in Fig0] except the fact that the degree of the suppression for semi-elliptic 
coupling spectral density is larger than the one for the Gaussian spectral density. 
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FIG. 5: Time evolution of I(t) for semi-elliptic coupling spectral density with r = 1, s 
7 P = 0.15; (a) without pulse application, (b) pulse interval t s = (c) pulse interval r s 
for pulse interval r s = t p . 

Defining the asymptotic peak value of I(t) for a pulse interval f s as P(f s ), we show 
the pulse-interval dependence of P(f s ) in FigJHJ The solid (dotted) line in FigJH| plots the 
P(f s ) for semi-elliptic (Gaussian) coupling spectral density. We find that P(f s ) has a local 
maximum when the pulse interval f a is close to 2/T. Since P(2tt) is nearly equal to P(l) as 
shown in FigEJ we find that the same degree of decoherence suppression is obtained for much 
longer pulse interval by paying attention to the dynamical motion of the reservoir. However, 
the effectiveness of the SPC decreases with increasing of the width 7 P of the coupling spectral 
density. 



= 3, and 

= T f: (d) 



B. Lorentzian coupling spectral density 

Now we assume the coupling spectral density to have the Lorentzian distribution as, 

h L (e) = \— 2 . (33) 
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FIG. 6: Pulse-interval dependence of the asymptotic peak value P (f s ) for the same parameters 
as in Fig. |I]and|2l The solid (dotted) line shows the P(f s ) for semi-elliptic (Gaussian) coupling 
spectral density. 

The Lorentzian coupling function has been often used in quantum optics. A relaxation 
process of an atomic system or quantum dots in a high-Q cavities has been described with 
a structured reservoir. The structure is determined by a distribution of coupling constants 

and often deseed with a Lo—functiooQQQ. 

We show the time evolution of I(t) for s = 3 and 7 P = 0.15 in Fig|7| Without 7r pulse 
application, I(t) shows the damped oscillation as in FigEfa). Contrary to the previous two 
cases of Gaussian and semi-elliptic distribution, FigE^b) shows that we cannot obtain the 
sufficient decoherence suppression for short pulse interval r s = Increasing the pulse 
interval to r s = y, we find that the degree of suppression becomes worse(FigJTfc)). The 
time dependence under SPC is shown in FigJT^d). We find that the almost the same time 
evolution as the one without pulse control. 

It should be noted that the ineffectiveness of the SPC does not come from the fact 
that the decay occurs faster for the Lorentzian coupling spectral density than the non- 
Lorentzian one. In order to make the point clear, we show the time dependence for smaller 
7 P (= 0.04), where we find the larger amplitude of oscillation in the time evolution without 
pulse application in FigJHJ The SPC makes the situation even worse, and is ineffective for the 
system with Lorentzian coupling spectral density. We will discuss the physical background 
of the ineffectiveness of the SPC in the next section. 
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FIG. 7: Time evolution of I(t) for Lorentzian coupling spectral density with j p = 0.15; (a) without 
pulse application, (b) pulse interval r s = 7^7, (c) pulse interval r s = (d) pulse interval t s = t p . 
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FIG. 8: Time evolution of 1(f) for Lorentzian coupling spectral density with j p = 0.04. The 
evolution with the SPC is compared the one without pulse application. 

IV. DISCUSSION 

Let us consider the reason why the SPC is ineffective for the case of Lorentzian coupling 
spectral density. Here we use a picture which provides us an qualitative understanding of the 
physical process of the SPC. The picture is obtained by replacing the original boson reservoir 
by a two-step structured reservoir where a single harmonic oscillator is coupled to a new 
"reservoir" with a coupling function which is different from the original one. (The "reservoir" 
means the newly introduced reservoir that causes the decay of the single harmonic oscillator.) 



The single harmonic oscillator is called as a quasi mode for atom-cavity system 



20 
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22, 
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251 ]. In the following, we call 



1231 . I24j or an interaction mode for electron-phonon system 
the new harmonic oscillator as the interaction mode(FigEJ). The motion of the interaction 
mode, which is determined by the coupling to the "reservoir", is characterized by the original 
coupling spectral density. The frequency (decay constant) of the motion of the interaction 
mode corresponds to the center frequency (width) of the original coupling spectral density, 
respectively. 




(b) 




FIG. 9: Schematic representation of the two pictures for the boson system: (a) the normal mode 
picture, (b)the interaction mode picture. 

The application of a tt pulse causes time reversal to the two-level system. Since the 
two-level system is coupled to the interaction mode, the degree of time reversal depends 
on the reversibility of the interaction mode. In the SPC, we equalize the pulse interval to 
the oscillation period of the interaction mode. If the reversibility of the interaction mode 
partially remains at the pulse application times, the SPC is effective for suppression of 
decoherence. 

When the original spin-boson interaction is characterized by the Lorentzian coupling spec- 
tral density, it has been shown that the interaction mode is coupled to the "reservoir" with 
a flat (white) co» P Hn g spectral de^Q Q. 

The interaction mode shows the Marko- 
vian nature, which indicates the irreversibility of the motion of the interaction mode (see 
Appendix). The SPC is ineffective for the Lorentzian spectral density. 

For the case of the non-Lorentzian coupling spectral density, we can also use the two-step 
structured reservoir where the interaction between the interaction mode and the "reser- 
voir" is characterized by a non-white coupling spectral density. This implies that the time 
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evolution of the interaction mode is non-Markovian, and partially reversible at the pulse 
application time. One should note that the non-Lorentzian coupling spectral density may 
not always guarantee the effectiveness of pulse control when the couplin g sp ectral density 
has a slow power-law tail at high frequencies as a recent study indicates [13]. However, in 
many systems such as localized-electron phonon systems, the shape of the coupling spectral 
density has clear maximum and minimum cut-off frequencies. The detailed characteristics 
of the spectral density for an effective pulse control deserve further research. 

We should remark that in the limit of short pulse interval, the pulse control is effec- 
tive even for the Lorenzian coupling function. In this case, the interaction mode does not 

nnn 

oscillate, and the result is consistent with the ordinary dynamical decoupling lid 111]. 



V. CONCLUDING REMARKS 

We have proposed a new strategy to suppress decoherence by multipulse control, which 
is done by synchronizing a 7r pulse train with the dynamical motion of reservoir. We have 
discussed the effectiveness of SPC for the non-Lorentzian and the Lorentzian coupling spec- 
tral density. For the former case, we find a periodic recovery of a quantum superposition at 
the pulse application times, whereas we cannot obtain the recovery for the latter case. 
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Our scheme is somewhat similar to the synchronized quantum-beat echoes(SQBE) 
in the sense that pulses are applied synchronously with the dynamical feature of reservoir. 
However, the SQBE is obtained by accumulating the response to each pulse whose area is 
much smaller than n. This means that the sufficient time reversal does not occur in the 
SQBE scheme and that the obtained echo is different from ours where essential physical 
origin is time reversibility caused by the each ir pulse application. 

The SPC provides us an another kind of method to suppress the decoherence by paying 
attention to the dynamical motion of the reservoir. We hope that the synchronized pulse 
control might extend the possibility of the pulse control of decoherence. Especially, the 
drastic alteration in the feature of the quantum superposition by changing the pulse interval 
might indicate that the pulse application plays both roles to keep and erase a quantum 
memory. 
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APPENDIX: INTRODUCTION OF THE TWO-STEP MODEL 

In this appendix, we briefly review how to introduce the two-step model with using the 
interaction mode(or the quasi mode). The interaction mode is the single harmonic oscillator, 
and the two-level system couples only to this oscillator. The annihilation operator of the 
interaction mode is defined by the linear combination of those for the original normal modes 
as 

B = g- lS ^h k e k b k , (A.l) 

k 

where 




(A.2) 



Then, Eq.flD can be rewritten as 

H SB = h\e)(e\ h Mbk + b\) = hg\e)(e\{B + (A.3) 

k 

The new "reservoir" modes are determined to be orthogonal to the interaction mode, and 
the many "reservoir" modes are orthogonal with each other. These boson systems other than 
the interaction mode do not couple to the spin system, and are called the "reservoir" modes. 
We should note that the "reservoir" here implies the environmental degrees of freedom for 
the interaction mode, and not for the spin system. The interaction mode and the "reservoir" 
are coupled to each other as 

n BR = h(B 9jR) + B j J2 9j R s)> ( A - 4 ) 

where Rj (i?j) is the annihilation (creation) operator of the oscillator for "reservoir" modes. 
The transformation from the original normal modes to the interaction and "reservoir" modes 
provides a new picture where a two-level system interacts with a single harmonic oscillator 
which contacts with a "reservoir". 
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When we consider a subsystem which consists of the interaction mode and the "reservoir" 
with T = 0, the Markovian master equation which describes time evolution of the interaction 
mode is solved to give the exponential decay of the coherent state amplitudes {26;], whose 
decay constant is determined by the coupling spectral density gj. The center frequency of 
the coupling spectral density in the original normal modes corresponds to the frequency of 
the interaction mode, and the width of the coupling spectral density is associated with the 
decay of the interaction mode arising from the coupling between the interaction mode and 
the reservoir. The Lorentzian coupling spectral density expressed by Eq. (I33J) implies a 
Markovian time evolution of the interaction mode, which indicates the irreversibility of the 
motion of the interaction mode. 

The quasi mode in the quantum optics has been obtained by an analogous procedure j^oi 
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22l |23j . We have another example to show that the two-step model with the white "reser 



voir" is equivalent to the normal mode picture with the Lorentzian coupling function 
Other than these examples, the two-step model has also been used to describe an effect 
of coupling between the nuclear reaction coordinate and the other coordinates on electron 
transfer in biomoleculesj^, Sj- 
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